Photo-excited Carrier Dynamics and Impact Excitation Cascade in Graphene 
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Photo-excitation in solids can trigger a cascade in which multiple particle-hole excitations are 
generated. We analyze the carrier multiplication cascade of impact excitation processes in graphene 
and show that the number of pair excitations has a strong dependence on doping, which makes 
carrier multiplication gate-tunable. We also predict that the number of excited pairs as well as the 
characteristic time of the cascade scale linearly with photo-excitation energy. These dependences, as 
well as sharply peaked angular distribution of pair excitations, provide clear experimental signatures 
of carrier multiplication. 



Converting light to electrical currents or voltages is 
a complex, multi-step process which involves photo- 
excited particle-hole pairs undergoing scattering by ambi- 
ent charge carriers, by other photoexcited carriers and by 
lattice vibrations. One of the key questions in the field of 
optoelectronics is identifying materials in which carrier 
multiplication can occur, i.e. a single absorbed photon 
yielding a large number of particle-hole pairs as a result 
of the primary photoexcited pair producing secondary 
pairs. Efficient carrier multiplication relies on a combina- 
tion of characteristics such as a wide band of states with 
a large phase space density for pair excitations, strong 
electron-electron scattering, and not too strong electron- 
phonon interaction. While graphene is by no means a 
unique example of a system with these properties, it is 
believed to fit the bill better than other materials. This 
has motivated an intense investigation of photoexcitation 
processes in graphene-based systems [TH12]. 

One aspect of graphene that distinguishes it from other 
materials is its truly two-dimensional structure which 
renders electronic states fully exposed. Photoexcitation 
in such a system generates photoexcited carriers that can 
in principle be extracted by a vertical transfer process, 
e.g. in a sandwich-type tunneling structure. Vertical car- 
rier extraction eliminates carrier loss in a lateral trans- 
port betwen photoexcitation region and contacts, often 
an important limiting factor for optoelectronic response 
in semiconductor systems. 

Despite intense interest, the photo-excitation cascade 
in graphene remains poorly understood. Theory predicts 
that the linear dispersion of charge carriers acquires a 
negative curvature due to electron-electron interactions, 
d 2 e(k)/dk 2 < [T3], which inhibits decay via electron- 
electron scattering in undoped graphene [11]. However, 
while the prediction of negative curvature appears to be 
in agreement with transport measurements [T3] , ARPES 
experiments support the notion of interaction-mediated 
decay [T51[TB]; interaction- induced quasiparticle decay re- 
mains the subject of ongoing debate [TT1 H21 IT7| . 

Besides its immediate utility for optolectronics, photo- 
excitation cascade in graphene is of interest because of 
its analogy with jets in particle physics, which are nar- 



(b) Branching Ratio, J/Jel-ph 




-12 -10 -8 -6 -4 -2 2 4 6 8 10 12 
Photoexcited Carrier Energy, t' jE F 

FIG. 1: (a) Impact excitation (IE) cascade of a photo-excited 
carrier with initial energy E®. Each cascade step involves 
electron-hole pair excitations with energy Ae ~ Ef, where 
Ef is Fermi energy. The net number of generated pairs and 
the relaxation rate depend strongly on Ef [see Eqs. 0,(16)], 
and thus can be tuned by gate voltage, (b) E nerg y-loss rate 
via impact excitation, J [blue curve, see Eq.( |16[ )], and op- 
tical phonon emission, J7el-ph HI] (red curve) for a typical 
doping of Ef = 0.2 eV. (Inset) Branching ratio, J/J c \- P h 
vs. e and Ef- (c) Transition probability (in units of h~ 1 ) 
obtained from Eq.(15l. IE processes with different initial 
energies, ei/E F = -12, -10, . . . , 10, 12 (here Ae = a - e'). 
Electron (hole) contributions shown by blue (red) curves. 



row cones of hardrons and other particles produced in 
high-energy particle detectors. Massless dispersion in 
graphene limits the kinematics of photoexcited carriers 
undergoing scattering by other carriers in much the same 
way as the kinematics of hadronization of quarks or glu- 
ons in a particle physics experiment, where jets are cre- 
ated. Because both the primary photoexcited carriers 
and the secondary carriers share the same linear disper- 
sion, collinear scattering is enhanced. As we will show, 
this results in sharp collimated angular features reminis- 
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cent of ultra-relativisitic jets in high-energy physics. 

Here we show that the photo-excitation cascade 
in doped graphene is distinct from that in undoped 
graphene. We identify impact excitation (IE) as the scat- 
tering process that dominates carrier relaxation dynam- 
ics in this system. Multiple secondary electron-hole (e-h) 
pairs produced by IE scattering involving a photo-excited 
carrier and ambient carriers in the Fermi-sea can lead 
to efficient carrier multiplication. Our analysis predicts 
that IE processes result in a chain-like cascade consisting 
of sequential steps with relatively small energy loss per 
step Ae ~ Ep, where Ef is the Fermi energy in graphene 
doped away from neutrality (see Figjlja) and (c)). As we 
shall see, both the number of pairs produced in the cas- 
cade (carrier multiplication factor) and the characteristic 
energy for the pairs are highly sensitive to doping. As a 
result, the key parameters of photo-excitation cascade in 
graphene are expected to be gate tunable in a wide range. 

As we shall see, the IE rate takes highest values al- 
lowed by unitarity, T ~ Ef/{2-kK). This fast char- 
acteristic rate makes this scattering process a highly 
efficient relaxation pathway which dominates over the 
phonon-mediated pathway in a wide range of energies 
(see Figflja). Indeed, typical dopings of Ep = 200 meV 
yield Ep/h = (2-7rx)0.48 x 10 14 Hz which corresponds to 
a scattering time of t w 20 fs (this is consistent with 
the values found for the transport mean free path due 
to carrier-carrier scattering extrapolated to high energies 
|19j). This is far faster than electron-optical phonon scat- 
tering rates [HI [19]. As a result, IE produces an energy 
relaxation rate, J , that dominates over the energy relax- 
ation rate from the emission of optical phonons, jTei-ph 
(see Fig. [ljb)) for typical dopings and photo-excited car- 
rier energies (see also Fig. [ljb) inset); IE controls the 
photo-excited carrier relaxation dynamics in graphene. 

The dependence on excitation energy Eq and Fermi 
energy Ep provides clear experimental signatures of this 
relaxation mechanism. In particular, the average number 
of e-h pairs produced in the cascade triggered by a single 
photo-excited electron, (N), is 



(N) = 



de 
(Ae) 



(Ae) = 



r(e) 



(1) 



where (Ae) is the average energy loss per step, J(e) 
and r(e) are the IE energy-relaxation and scattering 
rates respectively (see Eq.(16l). Here El ~ Ep is 



a low-energy cutoff corresponding to the energy below 
which IE processes are quenched; we use El = 2Ep 
(see discussion below). Fig{2^a) indicates that (N) ex- 
ceeds unity and grows quickly for Eq above few Ep (red 
curve). Since (Ae) ~ Ep, we find that (N) scales as 
Eo/Ep. In particular, an approximately linear depen- 
dence (N) « 0.55E Q /E F is found for E /E F > 1. 

Similarly, the time it takes for the photo-excited elec- 
tron to completely decay, At, also exhibits strong Eq and 
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FIG. 2: (a) Average net number of e-h pairs produced in the 
cascade triggered by a photo-excited electron with energy Eq 
(red curve). Cascade duration, At = G/(u [eV])fs, see Eq.|2| 
(blue curve), (b) Angular dependence for the e-h excitation 
rate T(ip), where ip is the angle between ki and q, see Fig[3|a). 



Ep dependence. This fairly short time, on the order of 
hundreds of femtoseconds, is 



At 



de G(Eo/E F ) 



J(e) E F [eV] 



fs. 



(2) 



where Q is a dimensionless scaling function (blue curve 
in Fig. [2ja)). As shown in Fig. [2](a), Q scales approx- 
imately linearly with Eo/Ep, yielding a At that scales 
linearly with the excitation energy. For a typical dop- 
ing value of Ep = 0.2 eV and initial photo-excited car- 
rier energy Eq — leV we find At « 0.12 ps, far faster 
than typical electron-lattice cooling time scales found in 
graphene 20 -22J. This separation of time scales means 
that the energy relaxation cascade occurs independently 
of electron-lattice cooling. 

Lastly, the angular distribution for impact excitation 
transitions is highly anisotropic. This produces a strong 
scarch-light-type structure peaked along the preferred di- 
rection of momentum transfer shown in Fig. [2] (b). 

Our system is described by the Hamiltonian for N = 4 
species of massless Dirac particles, 

H = <l J iA V<J ■ k )^M + ^l-el> ( 3 ) 
k,i 

n e i- e i = ± J2 ^(q)4 +q ,i4'- q , i V'k', j Vk,i.(4) 

q,k,k' 

Here i,j = 1...N and V(q) = 27re 2 /|q|K is the Coulomb 
interaction. Importantly, transitions in a massless Dirac 
band governed by the Hamiltonian ftty are subject to cer- 
tain kinematical constraints [TTJ H2J T These constraints 
arise due to the combined effect of linear dispersion in 
two Dirac cones, E±(p) = ±w|p|, and momentum con- 
serving character of carrier scattering. Here we ana- 
lyze the simplest case of a two-body collision. Each of 
the two particles participating in a collision can make 
transitions between states in the upper and lower Dirac 
cones which we denote by + and — . Two kinds of 
transitions can be distinguished: intra-band transitions 
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(H ► + or > — ) and inter-band transitions (H > — 

or ——>+). Since momentum change in any transi- 
tion satisfies ||pi| - |p 2 || < |pi - p 2 | < |pi| + |f>2 | , 
the intra-band transitions can only occur when the en- 
ergy and momentum change are related by |Ae| < v|Ap|, 
whereas the inter-band transitions are possible only when 
|Ae| > v\Ap\. 

The scattering process of interest, pictured in Figjj^a), 
involves a photo-excited carrier with high energy and mo- 
mentum e kl Ep, | kx | 3> kp, which is scattered to a 
lower energy state having momentum k' x with recoil mo- 
mentum q = ki — ki given to an electron in the Fermi 
sea. The latter process results in a particle-hole pair ex- 
citation, as depicted by a transition from k 2 to k 2 in Fig. 
[3] (a). The transition rate for this process, evaluated by 
the standard Golden Rule approach, takes the form 

W KM = ^ ]T / k2 (l-/ k ' 2 )F k2 , k ,^q| 2 (5) 

q, k 2,k 2 

X ^, k l+q^ k2 , k 2-q^( £ ki - £ kl + £lc 2 - £k 2 )- 

Here / k is a Fermi function, and F kk ' = |(k's'|ks)| 2 is 
the coherence factor (s, s' = ± label states in the electron 
and hole Dirac cones) . We treat the Coulomb interaction 
which mediates scattering between the photo-excited car- 
rier and the carriers in the Fermi sea by accounting for 
dynamical screening in the RPA approximation: 



v° 



£(w,q)' 



£( W ,q) = l-Vq°n(q, W ), 



(6) 



where V® = 27re 2 /|q|K and e(a;,q) describes dynamical 
screening. Here n is the polarization operator 

/(ek, s ) - /(ek+q, s ') 



k,s.s' 



w + e k , s — e k +q, S ' 



(7) 

with the band indices {s,s'} = ±. This includes both 
intra- (s = s') and inter- (s 7^ s') band contributions 

For Eq.([5]) to give a non- vanishing result, the tran- 
sitions ki — > ki, k 2 —> k 2 must occur in like pairs, 
both intra-band or both inter-band. Since the transi- 
tion ki — > ki is restricted to be within a single band, the 
transition k 2 — > k 2 must also be intra-band. As a result, 
the relaxation of the photo-excited carrier via inter-band 
scattering is blocked, whereas intra-band scattering is al- 
lowed. Kinematical blocking of inter-band processes can 
in principle be relieved by three-body (or, higher-order) 
collisions which are not discussed here. Such processes 
may become important at high excitation power, how- 
ever we expect the effect of such processes to be weak in 
the low excitation power regime discussed below. 

As shown below, the typical energy of an excited pair 
is much smaller than the photo-excitation energy e kl . 
Anticipating this result, it is convenient to factorize the 
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FIG. 3: (a) Kinematics of intraband carrier-carrier scattering 
in doped graphene. A photo-excited electron makes transi- 
tion from ki to ki by exciting an electron-hole pair from the 
Fermi sea from k 2 to k a . (b) Spectral function _R(q, oj) of 
particle-hole excitation as a function of momentum transfer 
and energy transfer per scattering event, (inset) Angular dis- 
tribution of normalized energy resolved transition rate (see 
for fixed values of oj/Ef = 0.2,0.5, 1, 10 (red, black, 
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Eq. _ 

green, and blue curves respectively) . The near-collinear char- 
acter of scattering at high cj is manifested in narrowing of the 
angular distribution. 



transition rate by expressing it through the spectrum 
of secondary pair excitations by following the standard 
procedure [23 . We write (5(6^ — e kl + e k ^ — e k2 ) — 
duj5{e-\ IL '_ L — e kl +w)r5(e k2 — e k2 —u). Next, we use the 
identity 

/k 2 (1 - /k 2 ) = (/k 2 - /k 2 ) x ( N(e K - e k2 ) + 1) , (8) 

where N(u) = l/(e w / fefsT — 1) is the Bose function taken 
at the electron temperature. Finally, we express the sum 
°f (/k 2 — /k 2 )^( e k 2 — £ k 2 — through a suitably defined 
susceptibility 

X"(q,w) - ^V^F kik+q (/ k -/ k+q )5(e k+q -e k -a;), (9) 



which can also be written as x"(q, w) = -jlmll(q,w). 
This yields a compact and intuitive expression for the 
total scattering rate: 

/oo 
dwP(«), (10) 

K i "°° 

P(w)=A^|V q | 2 F kliki x"(q,o;) ( 5(e kl -e kl +a;),(ll) 



where A = ^-[N{uj) + 1)][1 - /(e k - u)] and ki = ki q. 

As we show below, the typical energy and momentum 
transferred per scattering is of order Ep and Ep/v, re- 
spectively. These values are much smaller than those 
of the photo-excited electron. We can therefore approx- 
imate Ficj.i^ ~ 1, /(ffc — « and write the delta 
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function as 



function in energy to a delta function in angles in Eq.( 11 1, 



<5(£k' - £ki + w) ~ 5(«|q| cos V3 - u), 



(12) 



where is the angle between ki and q. The approx- 
imation |q| <§C |ki|, \uj\ -C w|ki| is appropriate under 
realistic conditions: for example, visible light frequencies 
translate to ej. = hf/2 — 750 meV which is considerably 
larger than Ep for typical doping values. Eq.(10l then 
yields the angle dependent transition rate 



T(ip) = / doj-f^ip), 



(13) 



d 2 q 

W) 2 



i?(q, w)<5(v|q| cos<^ — w), (14) 



where R(q_, u)) = A\Vq\ 2 x" (<l, . We evaluate the spec- 
tral function i?(q, u>) using the RPA-screened interaction, 
Eq.([6]),and susceptibility expressed through the polariza- 
tion function from Ref. [24] and the interaction param- 
eter a = e 2 /{nTiv) — 0.73. The angular distribution, 
r(<p), as well as the energy resolved distribution, 7^ (</?), 
feature interesting angular patterns (see Fig[2] (b) and 
Figj3] inset). Note in particular a sharp search-light-type 
structure corresponding the preferred direction of mo- 
mentum transfer, q, in the IE process. The peaks move 
closer to <p = as ui increases, indicating that the carrier- 
carrier scattering with a high energy transfer is nearly 
collinear. This is analogous to radiation pattern for an 
ultra-relativistic particle becoming focused along particle 
velocity [23], 

The angular patterns indicate that IE does not scatter 
into the angles \<p\ > ir/2 and vanishes at <p = 0. The ori- 
gin of these features can be understood from the spectral 
function i?(q, us) shown in Fig. [3] (b) which peaks close 
to ui = v\q\ but vanishes for both tp = arccos(a;/ti|q|) = 
and tp — arccos(w/v|q|) = ±7r/2. Scattering angles 
\<p\ > tt/2 correspond to negative ui, i.e. to the IE pro- 
cesses in which the photo-excited carrier gains energy 
from the Fermi sea. These processes are suppressed by 
the Bose factor in A. Naively one would expect that 
tp = 0, which corresponds to perfectly collinear scat- 
tering, would maximize transition rates since the phase 
space for particle-hole excitations in x" diverges when 
u|q| = u>. However, because the same divergence also ap- 
pears in II, the RPA-screened interaction makes i?(q, ui) 
vanish on the cone u> = v\q\. We note that this sup- 
pression of perfectly collinear scattering is highly sensi- 
tive to the approximation employed to describe dynam- 
ically screened interaction. The angular patterns F(<£>) 
and 7„ (tp) can therefore be used as a test of the validity 
of the RPA approximation. 

The same approach can be used to obtain the energy 
spectrum of pair excitations. In the following, however, 
we study the full energy dependence of P(Ae) not limit- 
ing ourselves to the asymptotic behavior at high photo- 
excited energies. Using a Jacobian to convert the delta 



we perform the angular integral in Eq.(ll I to obtain 



P(Ae) = f 
Jo 



2|ki| - (Ae/vh) - |q|cos^ 
(27r) 2 |ki||q|?iu simp 



R(q, Ae)qdq, 
(15) 

where tp is the angle between ki and q and satisfies 
(|k!| 2 - 2|k 1 ||q|cos ¥ > + Iql 2 ) 1 / 2 - |ki| = Ae/{vh). Nu- 
merically integrating Eq.(15l and taking e /c^T yields 
transition probabilities P(Ae) shown in Fig. [ljc) for dif- 
ferent initial photo-excited energies = . We find 
that P(Ae) peaks close to Ae sa Ep and decays rapidly 
for Ae ^> Ep. This non-monotonic dependence arises 
from the competition between the available phase space, 
which grows with Ae, and the Coulomb Interaction form 
factor, which decreases with |q|. 

The efficiency of IE scattering noted above can be 
linked to the large values of Ep in graphene. The relation 
between efficiency and Ep can be clarified by simple di- 
mensional analysis. We note that P(Ae) depends on Ae 
essentially via the dimensionless parameter x = Ae/Ep. 
This is clearly seen e.g. from pair excitation spectrum 
shown for different values of initial energy e^ = ei^ in 
FigjTJc): the width and profile of -P(Ae) has a very weak 
dependence on e^. This can be captured by writing the 



scattering rate V (Eq.(10)) as well as energy relaxation 
rate J = AeP(Ae)dAe in the form 



r(e) 



IT 



P(x)dx, 



e/E F 

xP(x)dx, 



(16) 



where we introduced dimensionless P(x) = KP(Ae). Us- 



ing -P(Ae) evaluated from Eq.(15), we obtain the fast 



impact excitation energy relaxation rate shown in Fig. 
[TJb) for Ep = 0.2 eV. Comparing J{e) from impact 
excitation with the energy relaxation rate arising from 
the emission of optical phonons, vTd-ph [18] . we obtain a 
branching ratio J /J e i-ph shown in the inset if Fig. ITVb) . 
While Jci-ph does not depend on carrier density [18], J 
does. As a result, gate voltage can be used to tune the 
branching ratio t7/j7" e i— ph by up to an order of magnitude 
as illustrated in the inset of Fi g. [l) (b) . 

The Ep dependences in Eq.(16l manifests in observ- 



ables such as the average number of secondary e-h pairs 
produced in a single photo-excitation cascade, (N) and 
its total cascade time, At. These quantities are related 
via (N) = L t Tdt. Using de/dt — —J{e) combined 
with Eq.(l6| we obtain Eq.Q for (N) and Eq.Q for 
the cascade time. In both cases, we have used a low en- 
ergy cutoff for the energy below which IE processes are 
quenched, El ~ Ep. Below the energy El, the relax- 
ation and scattering of the carrier from impact excitation 
slows dramatically and other relaxation processes dom- 
inate, for example energy relaxation via the emission of 
acoustic phonons. In evaluating Eqs.([T]),([2]), we used the 
value El = 2Ep below which the predicted value At 
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rapidly increases. The scaling of At and (N) with both 
excitation energy Eq and doping in Fig. [2]ja) are clear ex- 
perimental signatures of IE. Currently, particle-hole pair 
production and cascade times are the subject of intense 
experimental interest [2"tjH2"5] . 

In summary, the impact excitation scattering is pre- 
dicted to be a dominant carrier relaxation pathway in 
doped graphene. Multiple pair excitations are produced 
in a cascade triggered by a single photon, featuring an 
approximately linear scaling of the number of generated 
pairs and total cascade time with photo-excitation en- 
ergy. These dependences, as well as a sharply peaked an- 
gular distribution of e-h pairs, provide clear experimental 
signatures for IE-dominated cascade. Strong gate depen- 
dence of the cascade parameters affords a useful knob for 
the control of ultrafast scattering processes in graphene. 

We acknowledge financial support from the NSS pro- 
gram, Singapore (JS), the Office of Naval Research 
Grant No. N00014-09- 1-0724 (LL), and Fundacio Cellex 
Barcelona (KJT, FK). 
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APPENDIX A: ENERGY RELAXATION FROM 
THE EMISSION OF OPTICAL PHONONS 

An alternative channel for energy relaxation of photo- 
excited carriers occurs through the emission of optical 
phonons and gives an energy relaxation rate of J'ei-ph- 
The transition rate of this process [TO] can be described 
by Fermi's golden rule 

<; k ph = ^EiM^k)! 2 (17) 

q 

<5(Ae k ,, k + u; q )<5 k ,, k+q (iV(w q ) + 1), 

where Ae k / k = e k < — e k , w q = uiq = 200 meV is the 
optical phonon dispersion relation, and N(u)q) is a Bose 
function. Here k is the initial momentum of the photo- 
excited electron, k' is the momentum it gets scattered 
into, and q is the momentum of the optical phonon. The 
electron-phonon matrix element M(k',k) is 

|M(k',k)| 2 = 9o 2 F k , k ,, . 9o = -fi— (18) 

V2pw a 4 

where Fi^ic' is the coherence factor for graphene, go is 
the electron-optical phonon coupling constant |20j . p is 
graphene's mass density, and a is the distance between 
nearest neighbor carbon atoms. The energy- loss rate of 
the photo-excited carrier at energy e due to the emission 
of an optical phonon is 

Jel-ph(e) = Yl <7k Ph (4 - *) [1 - fM] . (19) 

k' 
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Integrating over q and k' we obtain 

Jci-ph(e) = %r-iv go[l-f(e-uio)](N(u>o) + l)v(e-Uo), 

(20) 

where v{e) — e/ '{2-kv 2 Ti 2 ) is the electron density of states 
in graphene. . Hence, l / e i-ph(e) varies linearly with the 
photo-excited carrier energy e > uj and vanishes for e < 
luq. Because the electron-phonon coupling with optical 
phonon is a constant, this result is to be expected from 
the increased phase space to scatter into at higher photo- 
excited carrier energy. 

To get an estimate of the energy relaxation rate, we 



estimate (N(u)q) + 1)«1 and 1 — /(e — oj ) w 0(e — ljq — 
E F ) and use p = 7.6 x 1CT 11 kgcm -2 and a = 1.42 A to 
obtain 

Jel-phCe) « e>i?F+ W o. (21) 

Using this we plot the energy relaxation rate for photo- 
excited electrons from the emission of optical phonons 
in Fig [l] (b). As shown in Fig. [l] (b), this channel is 
far smaller than the impact excitation channel described 
in the main text for typical dopings and photo-excited 
carrier energies. 



